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Non-commutative (NC) field theories can be mapped onto twisted matrix models. This mapping enables their 
Monte Carlo simulation, where the large N limit of the matrix models describes the continuum limit of NC field 
theory. First we present numeric results for 2d NC gauge theory of rank 1, which turns out to be renormalizable. 
The area law for the Wilson loop holds at small area, but at large area we observe a rotating phase, which 
corresponds to an Aharonov-Bohm effect. Next we investigate the NC 4 model in d = 3 and explore its phase 
diagram. Our results agree with a conjecture by Gubser and Sondhi in d = 4, who predicted that the ordered 
regime splits into a uniform phase and a phase dominated by stripe patterns. 



1. INTRODUCTION 

In the recent years there has been tremendous 
activity in the research of field theories on NC 
spaces. One field of application of such theories is 
the quantum Hall effect [Q. The boom of interest, 
however, was triggered by the observation that 
string and M theory at low energy correspond to 
NC field theories, see e.g. Ref. ||. 

The NC space is characterized by the non- 
commutativity of some of its coordinates. Here 
we are concerned with the case where only two 
coordinates do not commute, 

[x^x v ] = iOe^ . (1) 

For a review of field theory on NC spaces 
("NC field theory"), see Ref. f§. The non- 
commutativity yields non-locality in a range of 
O{y/0). This property implies conceptual prob- 
lems, but there is also hope that it might provide 
the crucial link to string theory or quantum grav- 
ity. From the technical point of view, it makes the 
perturbative renormalization even harder than in 
commutative field theory, because a new type of 
singularity occurs, which is UV and IR at the 
same time. Heuristically this can be understood 
immediately from relation ([!]) along with Heisen- 
berg's uncertainty relation. So far, there is no sys- 
tematic machine for absorbing this kind of mixed 
UV/IR divergences. 
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2. GAUGE THEORY ON A NC PLANE 

In particular, the NC rank 1 gauge action can 
be written in a form which looks similar to the 
familiar (commutative) case, 

S[A] = ^Jd d xT r (F^*F^) (2) 
= d^A v - d v A^ + ig{Afj, * A v - A v * A M ) , 
where the fields are multiplied by the star product 
f(x)*g(x) = e^^^L f(x)g(y)\ x=v . (3) 

This action is star gauge invariant. Also lattice 
formulations exist, but they can hardly be simu- 
lated: one would need star unitary link variables, 
construct a measure for them etc. 

In 1999 Ishibashi et al. observed g that NC 
U(n) gauge theories can be mapped onto spe- 
cific versions of the twisted Eguchi-Kawai model 
(TEK) Q , which has the action 

Stek[U] = -N/3 Z ^ ^U^IuUpl) , (4) 

Z^ = e 2mk ' L = Z* V)i {n<v,keZ,L = N 2 ' d ) . 

Since this system lives on one point there is only 
one N x N link variable in each direction. The 
factor Z^v is the twist. This mapping is based 
on Morita equivalence (an exact equivalence of 
the corresponding algebras) in the large N limit 
Qj, and on discrete Morita equivalence at finite 



2 



N, which corresponds to NC gauge theory on the 
lattice H . The equivalence holds if we choose 



N+l 



and it implies 9 



1 



-Na 2 



(5) 



where N must be odd, and a is the lattice spacing 
in the NC space. Also the Wilson loops, which 
are defined in the TEK in the obvious way (see 
below) are mapped onto the Wilson loops in NC 
gauge theory. 

This provides a non-perturbative access to NC 
gauge theory. The problems with the direct NC 
lattice formulation are all circumvented by this 
mapping onto a matrix model; hence this is a real 
counter example to the "no-free-lunch theorem" . 

Here we are interested in pure NC U(l) gauge 
theory in d = 2 and we summarize our numeric 
results of Ref. jjj. For analytical work on this 
model, see in particular Ref. || . 

Due to Eguchi-Kawai equivalence, the large N 
limit of the TEK at fixed j3 coincides with the pla- 
nar large N limit of commutative gauge theory, 
which has been solved analytically [||. In that 
case the Wilson loop obeys an exact area law, 



(W(IxJ)) = 



(6) 



1 



k(P) = -Mi-^) 



at 



From 
area" 



that framework we adapt the "physical 
a 2 IJ, where a = ^[k. However, for our 
purposes it is crucial to take the double scaling 
limit, which keeps Na 2 ~ N/f3 constant, so that 
the corresponding NC parameter 6 remains finite 
in the large N limit, c.f. eq. (||). 

This is different from the planar large N limit; 
the latter would lead to the Gross- Witten solu- 
tion, in agreement with the well-known equiva- 
lence of 8 — > oo and commutative gauge theory 
in the planar large N limit. 

The 2d TEK was simulated before but not 
with the parameters required here. Thanks to a 
trick introduced in Ref. Q we could apply a heat 
bath algorithm. 

We focus on the square-shaped Wilson loop, 

w 12 (ixi) = ^z^iuiuluyul 1 ) 



N 



(7) 



Note that W(I) := (W 12 (I x /)) is complex be- 
cause of the twist. In Fig. |l| we show its behavior 
in polar coordinates at N//3 — 32. We see that a 
double scaling limit does manifestly exist. Since it 
corresponds to the continuum limit of 2d NC rank 
1 gauge theory, we conclude that its Wilson loop 
is non-perturbatively renormalizable. At small 
area (with respect to 0(0)) it is practically real 
and follows the Gross- Witten area law of eq. (g). 
Beyond that regime |W(/)| does not decay any 
more, but the phase grows linearly in the area. 
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Figure 1. The polar coordinates of the complex 
Wilson loop W(I) plotted against the physical 



area A 



,2 t2 



At small areas it follows the 



Gross- Witten area law (solid line). At larger ar- 
eas the absolute value does not decay any more, 
but the phase increases linearly. 

Fig. H shows the asymptotic linear behavior of 
the phase <j) at different values of N//3 ~ tt9. In 



3 



the large area regime, we observe the simple re- 
lation 

= A/9 , (A = (al) 2 = physical area) (8) 

to a very high accuracy. We checked that 
it holds more generally for rectangular Wilson 
loops. Hence we can identify this law with the 
Aharonov-Bohm effect, if we formally introduce 
a constant magnetic field 



B = 1/9 



(9) 



across the plane. In fact, this corresponds exactly 
to the relation used by Seiberg and Witten when 
they mapped open strings in a constant gauge 
background onto NC gauge theory 0. Moreover, 
the very same relation was used in the description 
of the quantum Hall effect, where an electron in a 
layer was projected to the lowest Landau level [Q . 
Here we recover this law as a dynamical effect. 
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Figure 2. The phase of the Wilson loop W(I) 
at N = 125 for various values of N / (3, corre- 
sponding to different values of 9. At large areas 
A > 0(9) the phase <f> agrees with the formula 
4> = A/9, which is shown by straight lines. 

We fix again N/j3 = 32, and we consider now 
the connected Wilson 2-point function 



G 



(W) 



(W^(I x I)W vlt (I x I)) c e 



(10) 



(W) 



After a wave function renormalization G L 
j3^ a - 6 G'^ V ^ we observe also here a double scaling 
regime, see Fig. 3 (on top). 



Finally we consider the Polyakov lines 

P,(I) = Tr(C^) , P_ M (7) = Tr(£/t') . (11) 

Note that (P M (J)) = due to the phase symme- 
try (which corresponds to translation invariance), 
but the Polyakov multi-point functions are sensi- 
ble observables. Fig. 3 (below) shows that also 
the Polyakov 2-point function 



g[ p) {i) = (P M (/)P_ M (J)) 



(12) 



has a double scaling regime, if we apply the same 
wave function renormalization again, G. 
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Figure 3. The Wilson 2-point function G. 

top) and the Polyakov 2-point function G 2 (be 
low), at N = 195. After performing the same 
wave function renormalization we find a double 
scaling regime for each of the 2-point functions. 
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3. THE NC 4 MODEL IN d = 3 

We now summarize our ongoing work on the 
NC (j) A model in d = 3 We investigate 

the phase diagram and compare our observa- 
tions in particular with a conjecture by Gubser 
and Sondhi, which was obtained in d = 4 by a 
Hartree-type approximation They studied 

the phase-diagram in the m 2 / A 2 -6 A 2 plane for 
fixed coupling A, where A is a momentum cut-off. 
The conjecture can be summarized as follows: 

• For small 9 there is an Ising type phase tran- 
sition between a disordered and uniformly 
ordered phase. 

• For large enough 9 this transition changes 
its nature, due to UV/IR effects. This 
phase transition is driven by a mode of 
the scalar field with non-vanishing momen- 
tum, leading to a non-uniformly ordered or 
striped phase. This phase transition is first 
order in the regularized theory. 

• The theory is renormalizable within the 
one-loop self-consistent Hartree-type ap- 
proximation and the first order phase tran- 
sition turns into second order in the contin- 
uum limit. 



The action of the NC At/) 4 theory reads 



S 



d d x 



+ 



X 



*4 



(13) 



where the interaction term is a star-product of 
four scalar fields. We study this model in d = 
3, and we exclude the time direction from non- 
commutativity. The spatial coordinates satisfy 
eq. ©. 

A lattice version of this action can be con- 
structed, but as in Section 2 we map the sys- 
tem on a dimensionally reduced model to enable 
Monte-Carlo simulations. This mapping was al- 
ready carried out in Ref . Q . In our case the scalar 
field cj)(x, t) defined on a N 2 T lattice is mapped 
on N x N Hermitian matrices <fi(t). Their action 
takes the form 



S 



5 (r^(«)rt-£(t) 



(14) 



+i (j>(t + 1) - i{t))\ ^ 2 {t) + |0 4 (i) 



There are two kinetic terms according to the 
two underlying geometries. The matrices are 
called 'twist eaters' and act here as shift opera- 
tors. They are defined by the relation 



(15) 



where is the twist already introduced in Sec- 
tion 2. Again we have to use odd values of N. 

For this action we study the phase diagram, but 
in contrast to Gubser and Sondhi we consider the 
to 2 — A plane at fixed 9. By analogy we expect a 
striped phase for large enough A. 

Our results for N = 15 . . .45 are summarized 
in the phase diagram in Fig. |[ 
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Figure 4. The phase diagram of the 3d NC (p 4 
theory. The connected points show the separation 
line between disordered and ordered regime, and 
the hatched region marks the transition between 
uniformly ordered and striped phase at N = 35. 

So far we have been able to identify the separa- 
tion line between ordered regime and disordered 
phase and we observe a large N scaling. This 
separation is shown by the points connected with 
lines. The vertical lines indicate the transition 
region between uniformly ordered and striped 
phase at N — 35. For N < 35 we find the tran- 
sition between the two ordered regimes at larger 
values of A^A. Below we discuss the procedure 
that led us to this phase diagram. 
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We use the momentum dependent order 
parameter (M(k)), where 



iV J p — K ' 



(16) 



and 4>(p, i) is the spatial Fourier transform of 
(j>(x,t) ^ In particular (M(0)) reduces to the 
standard order parameter for the spontaneous 
breakdown of a Z<i symmetry. To find the separa- 
tion line between ordered and disordered phase, 
we decreased m 2 at constant values of A, starting 
in the disordered phase. Two examples at A = 35 
are shown in Fig. 
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Figure 5. XTie momentum dependent order pa- 
rameter against Nm 2 . Above we fixed NX = 2, 
which leads to a uniformly ordered phase, and be- 
low at NX — 20 leading to a striped phase. 

In these plots we show the standard order pa- 
rameter (M(0)) and the staggered order param- 
eter (M(l)). In the upper plot of Fig. || only 
(M (0)) is non-trivial below some value of m 2 < 0. 

2 Note that <f>(p, t) coincides with the Fourier transform in 
the commutative case. 



Therefore we are in the uniformly ordered phase. 
In the lower plot of Fig. ||, (M(0)) practically 
vanishes for all m 2 , but (M(l)) is non-zero, indi- 
cating a striped phase. Such measurements allow 
us to identify the separation region between uni- 
formly ordered and striped phase. 

To localize the separation line between ordered 
and disordered phase we computed the connected 
part of the two-point function of M(k). This 
function has a peak at the phase transition. Fig. 
H illustrates some examples for these measure- 
ments. 
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Figure 6. The two-point function of the quantity 
M{k), defined in eq. (f/6|j ; against Nm 2 at NX = 
2 (above) and at NX — 20 (below). 

In Fig. ^ we show some snapshots of configura- 
tions 4>(x, t) at N = 35. The dotted areas mark 
the regions where <f> > 0, and the blank areas rep- 
resent 4> < 0. At A A = 2 and for small |Am 2 | 
(Fig. 0a) we find such areas scattered all over the 
volume, and at Nm 2 -C (Fig. f?]b) (f>.(x, t) is ei- 
ther positive or negative for all x, as one expects 
for an Ising type phase. At AA = 20 we also find a 
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disordered phase (Fig. 0c), but for Nm 2 <C the 
uniformly ordered phase is replaced by a striped 
phase (Fig. 0d). 
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Figure 7. Snapshots of single configurations 
4>{x,t) at a certain time t, for N = 35. 

4. CONCLUSIONS 

We presented numeric results for NC field the- 
ory, which could be obtained through the map- 
ping onto twisted matrix models. 

In the first part we summarized the results of 
Ref. Q. Specific types of the TEK correspond 
to 2d NC U{1) gauge theory on the lattice. We 
approached the large N limit so that the non- 
commutativity parameter 9 is kept constant (dou- 
ble scaling limit). This describes the continuum 
limit of NC gauge theory. We found finite lim- 
its for the 1- and 2-point Wilson loop, as well 
as the 2-point Polyakov line, which confirms that 
this model is non-perturbatively renormalizable. 
This is the first evidence for the non-perturbative 
renormalizability of a NC field theory. 

At small areas, the Wilson loop follows an 
area law as in the (commutative) planar theory. 
At larger areas it deviates and reveals a new 



double scaling limit. This is a manifestation of 
non-perturbative UV/IR mixing. In that regime, 
the phase of the Wilson loop is given simply by 
the area divided by 9. This corresponds to an 
Aharonov-Bohm effect if we identify 9 with a 
magnetic field B = 1/9, as it is also done in the 
solid state jjjj and string || literature. 

We further studied the phase diagram of a 3d 
NC <fi 4 theory 12 . The separation line between 
ordered and disordered phase is localized and we 
observe large N scaling. In agreement with the 
qualitative conjecture in Ref. ||L3|] , the ordered 
regime splits into a uniform and a striped phase. 
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